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The influence of the salts KC1, NaCl, and Nal at molar concentrations on the a-helical folding 
kinetics of the alanine-based oligopeptide Ace-AEAAAKEAAAKA-Nme is investigated by means of 
(explicit-water) molecular dynamics simulations and a diffusional analysis. The mean first passage 
times for folding and unfolding are found to be highly salt-specific. In particular, the folding times 
increase about one order of magnitude for the sodium salts. The drastic slowing down can be traced 
back to long-lived, compact configurations of the partially folded peptide, in which sodium ions are 
tightly bound by several carbonyl and carboxylate groups. This multiple trapping is found to lead to 
a non-exponential residence time distribution of the cations in the first solvation shell of the peptide, 
j /^ ' The analysis of a-helical folding in the framework of diffusion in a reduced (one-dimensional) free 

energy landscape further shows that the salt not only specifically modifies equilibrium properties, 
but also induces kinetic barriers due to individual ion binding. In the sodium salts, for instance, 
the peptide's configurational mobility (or "diffusivity") can decrease about one order of magnitude. 
This study demonstrates the highly specific action of ions and highlights the intimate coupling of 
intramolecular friction and solvent effects in protein folding. 
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I. INTRODUCTION 



The complex problem of protein folding is typically interpreted in terms of a diffusive search through an effective, 
low-dimensional free energy landscape, where most of the countless degrees of freedom of the system have been 
integrated outii~— In this picture the effective diffusivity (or friction) arising from intrapeptide and peptide-solvent 
interactions as well as from orthogonal degrees of freedom varies considerably along the reaction coordinate (RC)ii£ 
In general, the internal friction constitutes a major contribution to the friction, such that solvent viscous drag and 
solvent-peptide interactions are not the only mechanisms that govern protein kinetics. 9,10 Internal and solvent-induced 
— , friction processes may even be intimately connected as was demonstrated for the loop formation rate of unfolded 
peptides, where the strongly denaturing salt guanidine hydrochloride modifies internal friction by specific binding 
mechanisms i 1 1 1 1 2 

The investigation of the salt-specific action on proteins in general ('Hofmeister effects') has a long history^ but 
the underlying mechanisms are still under exploration. In a recent series of experiments, for example, it has been 
shown that even the simple cations sodium (Na + ) and potassium (K + ) exhibit considerably different behavior in the 
interaction with protein surfaces, where Na + is favored over K + ji2r— One apparent reason is the stronger attraction 
of sodium to acidic (anionic) surface groups, in particular to carbonyl and side chain carboxylate groups. While these 
static properties have received much attention lately^! not much is known about their consequences to biomolecular 
kinetics. Experimental hints have been given in studies of Na + - and K + -specific polyglutamic acid aggregation 
kinetics^ folding kinetics of halophilic ('salt-loving' and very acidic) protein a 19 i 20 or DNA^i However, a detailed 
molecular understanding of the ion-specific action on biomolecular folding and assembly kinetics is still lacking. 

In particular for large concentrations (c > 1 M) salt effects are pronounced, highly sequence and salt-type specific, 
and typically lead to changes in protein solubility, stability, and/or denaturation that result in the so-called Hofmeister 
series for the precipitation of proteinsJ^ Although an order of magnitude higher than at typical physiological conditions 
(c ~ 0.2 M), large salt concentrations play a central biochemical role in the broad field of protein crystallization^ 
in food industry as fermentation additives^ and for the function and stability of biotechnologically interesting 
halophilic (salt-loving) enzymes^ Additionally, the study of protein structures in solvents with high salinity solvents 
is instructive as salt-specific effects are amplified and, important from a computational perspective, can be sampled 
more efficiently in molecular dynamics (MD) simulations, a popular tool nowadays for the theoretical study of protein 
folding, function, and stability^ 

An ubiquitous and fundamental secondary structure element of proteins is the a-helix, which is stabilized by 
(i, i + backbone hydrogen bonds involving four amino acids (aa) per turn. The majority of short (< 20 aa) isolated 
helices derived from proteins are unstable in solution, unless specific side-chain interactions stabilize them. Among 
those it has been demonstrated that alanine-based peptides have the strongest intrinsic helix propensity^ - — A very 
instructive model peptide was experimentally introduced by Marqusee and Baldwin 26 , showed that alanine-based 
oligopeptides with an A(EAAAK)„A pattern display a high cv-helix propensity, which is probably stabilized by Glu~ 
(E) and Lys + (K) salt bridges along the folded peptide. Marqusee and Baldwin also found that the simple salt NaCl 
has a destabilizing effect on the a-helical configurations of these peptides. 

Indeed, in a recent explicit-water molecular dynamics (MD) computer simulation^ the structural behavior of the 
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oligopeptide Ace-AEAAAKEAAAKA-Nme, in the following named the 'EK' peptide, has been explored in detail, and 
the stabilizing and destabilizing mechanisms in various highly concentrated aqueous electrolyte solutions have been 
identified. Among those it has been found that sodium (Na + ) destabilizes the helical structure stronger than potassium 
(K + ); it was also previously recognized that the reason for the destabilization of the salt bridges and the secondary- 
structure forming hydrogen lies in the higher affinity of Na + to carboxylate and carbonyl groups. In a recent follow-up 
paper the folding kinetics of the EK-peptide in pure aqueous solution was investigated and interpreted in terms of 
diffusion in a reduced (one-dimensional) free energy landscape involving a local coordinate-dependent diffusivity^ 

The aim of this work is to extend these previous studies on the 'EK'-peptide to the investigation of the specific 
effects of the salts KC1, NaCl, and Nal on peptide a-helical folding and unfolding kinetics. In particular, we explore 
how salts alter the mean folding and unfolding times and look for possible molecular reasons. Strongly salt-specific 
effects are found in the mean folding times, which can be attributed to the binding of individual ions to multiple, 
anionic peptide groups inducing transient cross-links between peptide fragments. In consequence, not only equilibrium 
distributions of configurations are modified, but also new, slow time scales in the peptide's configurational mobility 
arise due to enhanced internal friction. Salt effects arc thus reflected in both modified free energy landscape and local 
changes of the effective diffusivity. Our study demonstrates the highly specific action on protein folding kinetics by 
the individual binding of ions, and, more generally, exemplifies the intimate coupling between solvent and intrapeptide 
friction effects in protein folding. We believe that these mechanisms could be of general importance and transferable 
to a variety of biomolecular and polyelectrolyte systems. 

II. METHODS 
A. MD simulations 

Our MD simulations are performed using the parallel module sander.MPI in the simulation package Amber9.0 with 
the ff03 force field for the peptides and the rigid and nonpolarizable TIP3P water for the solvent^ All simulated 
systems are maintained at a fixed pressure of P = 1 bar and a temperature T = 300 K by coupling to a Berendsen 
barostat and Langevin thermostat^ respectively. The cubic and periodically repeated simulation box of edge length 
L ~ 36 A includes approximately 1500 water molecules. Electrostatic interactions are calculated by particle mesh 
Ewald summation and all real-space interactions (electrostatic and van der Waals) have a cut-off of 9 A. The peptide 
is generated using the tleap tool in the Amber packaged 

We investigate the helical folding and unfolding behavior of a 12 amino acids long peptide with the acetyl (Ace) 
and amine (Nme) capped sequence Ace-AEAAAKEAAAKA-Nme. This peptide can form three a-helical turns in the 
fully folded state, where Glu2 and Lys6, and Glu7 and Lysll, are potentially able to form a salt bridge, respectively^ 
The influence on a-helix folding kinetics of a large concentration of ~ 3.6±0.1 M of the simple monovalent salts NaCl, 
KC1, and Nal is investigated. We have simulated the system without salt for ~ 1.35 fis and including salt for ~ 2 fis 
for each salt type. We note here that the free energy along the RC q (see below) derived from replica-exchange MD 
trajectories at T — 300 K for the salt-free system gives very good agreement to the brute force approach^ indicating 
a decent statistical sampling by our trajectories. 

The considered salt concentrations result from 90 ion pairs in the simulation box. Cations and anions are modeled 
as nonpolarizable Lennard-Jones spheres with charge and interaction parameters as supplied by Dang^ as the default 
Amber parameters are known to be faulty»2i The Dang parameters show reasonable bulk thermodynamic properties 
in SPC/E water even for high concentrations^ Comparative calculations (unpublished) in TIP3P water show only 
small differences in hydration structure, and no qualitative difference in the binding to peptide groups^. The 
parameters used are summarized in previous work on the equilibrium structure of the EK peptide in salt at a different 
temperature^ We are aware of the weaknesses of ionic MD force fields for quantitative predictions in biomolecular 
simulations; this issue has already been discussed in our previous study, where however a reasonable description of 
helicity and destabilization with NaCl was observed when compared to experiments. Since the destabilization seems 
to be overemphasized for NaCl, we do not claim to be quantitative in our work, but focus on the discussion of the 
main effects and qualitative trends with the addition of salt; we believe these to be insensitive to the particular force 
field and relevant for a variety of experimental observations. 

B. Helicity and reaction coordinate 

Trajectory analysis is performed using the ptraj tool in the Amber9.0 packaged The helicity, i.e., the a-helical 
fraction, is identified using the DSSP method by Kabsch and Sander >2Z We focus on one RC Q, which is defined as the 
root mean square distance from a fully helical reference structure (i.e., with helicity equals to one), averaged over all 
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atoms of the peptide and which thus measures the deviation from the 'native' state; this quantity was previously found 
to act as an adequate dynamic RC in the salt-free case£ We note that different choices of the reference structure were 
tested and resulted only in small differences in the energy landscapes, i.e., local variations on the order of fractions 
of the thermal energy ksT; see the Supporting Information (SI) for additional details. Similarly, the RC-trajectories 
only marginally change when varying the reference structure and therefore yield the same kinetic behavior; examples 
of such trajectories are also found in the SI. Trajectories are recorded with a resolution of 20 ps, giving a total of ~ 
67500 data points for the simulations without salt and roughly 100000 data points for the runs including salt. For 
convenience, we define a rescaled RC by q — (Q — Q min )/(Q maK ~ Q mm ) such that the minimal and maximal values 
of the data points, denoted as Q min and Q max , are projected on the RC values q = and q = 1, respectively. The 
absolute minimum and maximum RC values are similar for all systems and are Q min ~ 1.0 A and Q max ~ 8.0 A. 



C. Bulk shear viscosities 



In order to get a more complete picture of the solvent properties we calculate bulk shear viscosities for NaCl, KC1, 
and Nal at the relevant concentrations in TIP3P water. We employ the same ionic force fields as above, but perform 
the simulations with the GROMACS 4.0 38 i 39 package due to the implemented viscosity calculation methods. In these 
simulations, the periodically repeated box has an edge length of L ~ 4 nm, with a total number of A w = 1910 
water molecules and A; = 135 ion pairs. For the pure water simulation we use A w = 2180 water molecules. After 
TVPT-equilibration we proceed with an A^T-production run of 50 ns. We compare two approaches to calculate the 
viscosity: firstly, we employ the Green-Kubo (GK) formul a 40 ' 41 

V f°° 

11 = foT J Q ( P -(*") P -( < o + i ))*o d *. (!) 

involving the off-diagonal components of the pressure tensor only. We obtained averaged viscosities over a correlation 
time of 5 to 20 ps. The latter choice reflects the fact that the viscosity converges rapidly, but exhibits large statistical 
errors for long correlation times.— 

Secondly, we perform a non-equilibrium perturbation method.— In this scheme, an external force is applied in the 
Al^T-simulation with the periodic acceleration profile 

a x {z) — Acos(fcz), (2) 

with k = L being the edge length of the box. The amplitude A should be chosen small enough in order not to 
drive the system out of the linear response regime and at the same time large enough to get good statistics. For a 
more detailed discussion we refer to previous work^ - and set the amplitude to A — 0.02 nm ps~ 2 . We then obtain the 
viscosity by calculating the average velocity profile of all particles^ - 

With the GK formula we find values of r/o — 0.31 ±0.01 10~ 3 kg/(m-s) for pure TIP3P water, corroborating 
with previous studiesi 3 - and r) = 0.58, 0.74, and 0.6±0.01 10~ 3 kg/(m-s)for the KC1, NaCl, and Nal solutions at 
a concentration of 3.6 mol/1, respectively. The periodic perturbation method yields the same results within a 5% 
error range. Compared to experimental values^ 4 - the MD simulation considerably overemphasizes the increase of 
the viscosity at this elevated salt concentration; indeed, the viscosity was experimentally found to increase by only 
roughly 5% for KC1 and 30-40 % for NaCl and Nal compared to pure water. This failure in describing the correct 
bulk viscosities of the electrolyte solutions must be attributed to inaccuracies in the force field. Note that the value 
for pure TIP3P water already considerably deviates from the experimental value (0.893 -10 -3 kg/(m-s) at 298.15 K) 45 
by more than a factor of two. 



D. Free energy and diffusivity profiles 

We assume that the stochastic time evolution of the probability ^(q,t) of finding a configuration with RC value q 
at time t is described by the one-dimensional Fokker-Planck (FP) equation 4 ^ 

= ^D(q)e-^^(q,t)e^\ (3) 

where D{q) is the (in general g-dependent) diffusivity, (3 = l/(k-QT), and F(q) = —kBT\n(^(q)) is the free energy 
obtained by Boltzmann-inversion of the time-averaged probability distribution (ty(q)). We employ a recently in- 
troduced method to estimate the diffusivity profile that takes advantage of the relation between average transition 
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times between different points along the RC, the free energy landscape F(q), and the diffusivity profile D(q). 8 The 
round-trip time defined as 

TRT(q, <?t) = sign(g - q t )[Ti p (q, q t ) + T ip (q t , q)}, (4) 

specifies the average time needed for an excursion starting at q, reaching q t at least once, and returning to q for 
the first time; mean first passage times for transitions starting at q and reaching qt are denoted by Tf v (q, qt). For a 
diffusive process described by the FP equation ([3]) the round-trip time is given by the integral 

f 9 i 1 

TRT(q,qt)= dq , (5) 
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which can be derived from the expressions for the mean first passage time given before^. Though Eq. ([5]) can in 
principle be inverted to obtain the diffusivity profile D(q) from the slope of the round-trip time curves, we choose 
a complementary analysis method within the present work to avoid artifacts due to insufficient statistical sampling. 
The FP approach assumes an underlying Markovian process and — as is easily seen from Eq. ([5]) — round-trip time 
curves (as a function of q) for different target points q t therefore only differ by a constant 

trt(<2, q t ) = r RT ((7, q[) + r RT (^, q t ). (6) 

The assumption of Markovian behavior however generally breaks down at short times and for unsuitable RCs, i.e., 
RCs that do not single out the transition state; see also previous literature on that topic— 

In our analysis of the simulation time series q(t) we discretize the RC in N = 50 intervals centered around q^ = 
(2i— 1)/100, i G {1, 2, . . . , 50}; mean first passage times between all possible pairs of bins are extracted from simulation 
data and converted into round-trip times using Eq. Q. To simplify the analysis, we assume a flat diffusivity within 
each of the following regions: (i) values of the RC q < <? 2 /3 corresponding to an almost perfectly folded helix, (ii) one 
or two unfolded helical turns equivalent to ^2/3 _! Q < <7i/3j and (iii) mostly unfolded states characterized by RC values 
q > qi/3: where the indices 1/3 and 2/3 denote the average helicity at these q- values. The values of the diffusivity in 
those three regions are used as fit parameters in Eq. ([5]) to best reproduce the round-trip time curves obtained from 
simulation data; the integral in Eq. ([5]) is computed numerically by linear interpolation of {('4 , (gW))}|£ 1 . Best fits to 
the round-trip data for different target points qt allow to determine the diffusivity (including an error estimate) for 
each of the three regions. Alternatively, the diffusivity profile can also be obtained by fitting to the average round-trip 
time 
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f RT (gW)E-^r RT (gW,g»), (7) 



which is less affected by statistical noise and which according to Eq. ([5]) is just shifted vertically w.r.t. the round-trip 
curve trt((7, <7t) for a specific target position q t . The results of both fitting procedures are shown and discussed in 
Sec. 1111 CI additional information concerning the fitting is provided in the SI. 



III. RESULTS AND DISCUSSION 
A. Equilibrium free energy landscapes 

The RC time series q(t) of all investigated systems is shown in Fig. 1: numerous folding- unfolding transitions 
on different time scales are discernible. Already at a first glance longer time scales (and fewer unfolding-folding 
transitions) seem to be involved in the systems with the sodium salts, which will be quantified and discussed in more 
detail later. Free energy profiles F(q) of all investigated systems are extracted from the trajectories q(t) and plotted 
in Fig. 2 together with the (time-averaged) helicity resolved along the RC q. In the salt-free case (black curve) in 
Fig. 2 (a), three minima are clearly visible: one at qi ~ 0.1, a second at q^ ~ 0.32, and a third, shallow one at 
<73 = 0.58. Judging from the helicity vs. q, plotted in Fig. 2 (b), these minima seem to correspond mainly to 1) a full 
helix with three turns, 2) a partially unfolded helix with two neighboring turns, and 3) one helical turn as well as fully 
unfolded states. Representative simulation snapshots taken at corresponding ^-values (i = 1,2,3) confirm this view 
and are also shown in Fig. 2 (b). We note that the value of which in absolute units corresponds to Qi ~ 1.6 A, 
deviates from (the reference state) due to thermal fluctuations. 

For KC1, the free energy F(q) is slightly shifted to favor unfolded states at larger q- values (q^ ~ 0.67), while the 
main features remain the same as without salt. The total helicity decreases slightly from 62% to 55%. For the sodium 
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salts, NaCl and Nal, the fully folded state at q\ becomes metastable, and the distribution is strongly shifted to the 
partially and fully unfolded states. In particular, the third minimum at q% deepens and broadens. Also in the systems 
with salt the minima mostly match with partially folded a-helical states as can be judged from the helicity vs. q 
curve. All total helicities and positions of the minima in F(q) are summarized in Tab. I. 

In the NaCl and Nal salts the total helicity significantly decreases to about 39% and 34%, respectively. The main 
causes have been discovered before^ 1 , and are twofold: firstly, specific Na + -binding to the glutamic acid side chain 
carboxylates interferes with E _ -K + salt bridging, and secondly, the binding of cations to backbone carbonyls perturbs 
secondary-structure forming hydrogen bondsi^I In comparison to Na + , the specific binding of K + is weak and KC1 
thus behaves less destabilizing. We note that quantitatively these effects are force field dependent while the binding 
trends have been corroborated in many experiments and various simulation studies^ - — For completeness, we plot 
radial distribution functions (rdfs) between the cations and the oxygen atoms from the backbone carbonyls and side 
chain carboxylates in Fig. 3. The dominance of sodium over potassium binding to the anionic peptide groups is clearly 
observable. 

TABLE I: Total helicity of the investigated systems and distinct positions in the free energy landscape shown in Fig. 2: locations 
of local minima are denoted by qi, i = 1, 2, 3, and values of the RC, where the average helicity is 2/3 or 1/3, by 92/3 and 51/3 
respectively. 



system 


helicity 




<72 


<73 


92/3 


<?l/3 


no salt 


0.62 


0.11 


0.31 


0.57 


0.34 


0.56 


KC1 


0.55 


0.10 


0.32 


0.67 


0.35 


0.65 


NaCl 


0.39 


0.08 


0.28 


0.59 


0.30 


0.52 


Nal 


0.34 


0.09 


0.31 


0.65 


0.32 


0.53 



B. Long-lived structures and specific ion binding 

Previous work on the EK-peptide stability in salt solution^! revealed that due to Na + binding long-lived loop- 
forming configurations occupy the region of intermediate to large q- values, in addition to one-turn and random coil 
states also present without salt. In these looped configurations a single Na + -ion is collectively trapped by a few peptide 
backbone carbonyls and side chain carboxylates. This leads to a partial wrapping of the oligomeric backbone around 
the ion. From a superficial inspection of our trajectories we find that these long-lived configurations involving tightly 
bound Na + -ions can be stable on a ~ 1-10 ns time scale. Representative simulation snapshots which are dynamically 
selected such that they exist for more than 2 ns are shown in Fig. 4: in panel (a) a sodium ion is bound and wrapped 
by the central part of the peptide, while the terminal parts exhibit partial helical turns. The bound sodium 'locks' 
this structure on times up to several nanoseconds. Similar long-lived situations are found, where two neighboring 
turns are correctly folded, while the rest of the peptide forms a loop around the cation, see Fig. 4 (b). Often also the 
glutamic acid side chain is involved as displayed in Figs. 4 (c) and (d) where the cation not only binds to a backbone 
carbonyl but also to the headgroup of the E~ side chain. In all these situations the peptide is relatively compact with 
one or two turns correctly folded such that the relevant g-region of these states is between 0.25 < q < 0.65, where 
the helicity is mostly between 2/3 and 1/3 (cf. Fig. 2 (b)). Configurations of this type have not or only rarely been 
found involving anions or a K + -ion, so that their existence must be attributed to the relatively strong binding of the 
Na + ion. 

The rigorous life-time analysis of the occurring configurations involving trapped ions is difficult to perform due 
to the variety and complexity of the somehow amorphous structures. We therefore resort to the analysis of cation 
binding times. In the inset to Fig. 3, we plot the binding time distribution Pb(t) of the cations K + and Na + in the 
first solvation shell of carboxylates and carbonyls (defined by the location of the first minimum in the cation rdfs in 
Fig. 3). While the distribution decays exponentially for K + with a time constant of about 50 ps, we find a much 
slower, nonexponential behavior for Na + which can be best fitted by a stretched exponential of the form oc exp(— t/r)P 
with t ~ 50 ps and j3 — 0.55. This indicates that long binding times on the order of nanoseconds are indeed possible 
and corroborate the existence of 'trapped' unfolded configurations in which peptide parts tightly wrap around the 
cation. 

We note that it is indeed well established that systems with multiple trapping or other manifestations of disorder can 
lead to anomalous kinetics. 47 In our investigated systems, the appearance and magnitude of trapping is controlled by 
the nature of the ions. The observed 'stretched' exponentials resemble the slow relaxation in glass-forming liquids. 48 



6 



Anomalous kinetics in peptide and protein dynamics have been indeed observed in simulations and experiments for 
certain dynamic variables or RCs^2r— and obviously question the general validity of diffusive approaches to predict 
long time dynamics in protein folding. While this complex issue is still awaiting resolution we proceed in this work 
with the interpretation of helix folding in the framework of simple diffusion; since the average folding / unfolding 
times (> 10 ns) are typically much larger then the ion binding times (~ 50 ps), we expect the long-term dynamics to 
be adequately described in terms of memory less diffusion in a free energy landscape. 

C. Folding kinetics and diffusional analysis 

Let us first analyze mean folding and unfolding times. In Fig. 5. (a) we plot the mean first passage time for folding 
from q > qi to qi given by Tf(q,qi) = T{ p (q, gi); the salt specific values of qi are found in Tab. U Without salt the 
typical folding time is about n ~ 20 — 30 ns in the region q > 0.3 before it quickly drops down to for q- values closely 
approaching q\. In KC1 Tf(g, qi) increases by a factor of about 2, while the sodium salts lead to a considerable slowing 
down of folding by one order of magnitude. The unfolding times T u f(g, g 3 ) (q < 53), cf. Fig. 5 (b), show a bit less 
variation between the salts. Without salt the typical unfolding time is about 30-40 ns while it may rise by a factor 
of 2-3 in NaCl or KC1. Note that for KC1 the unfolding is considerably slower than in the salt-free case although the 
free energy landscape is very similar. 

To get a grasp on the folding kinetics involving fewer helical turns we have also analyzed Tf(g, q2) for q > qi (folding 
by 1 or 2 turns to the 2-turn state) and r u f(g, (72) for q < q2 (unfolding by 1 turn to the 2-turn state), which are 
shown in Fig. 6 (a) and (b), respectively. The folding times decrease by a factor of about 2 when compared to the 
q — > qi folding, while the trends with salt remain the same, in particular a one order of magnitude slower folding in 
the sodium salts. Unfolding times T u f (g, 52) are relatively small and found between 3-7 ns. The variation between the 
salts is again less pronounced for unfolding when compared to the folding times and show a different ordering. 

Note that the mean first passage times shown in Figs. 5 and 6 arc in parts subject to substantial noise due to 
insufficient statistical sampling: for example, clear deviations from a monotonously increasing function (expected 
for diffusive dynamics) are observed in the KCl-folding times shown in Fig. 5 (a); this irregularity has its origin in 
an exceptionally long-lived (relative to the full trajectory duration) state of a specific peptide configuration in the 
KCl-time series in Fig. 1 within the time range 200 ns < t < 400 ns. Without going too much into detail, this specific 
configuration is characterized by two helical turns and the Ace-cap being buried between the hydrophobic side chains 
of one Lys and one Ala side chain. This configuration (not involving bound ions) is found rather frequently in all 
trajectories but with typically much shorter life times. 

Let us now turn to the interpretation of the mean folding times by the free energy F(q) and diffusivity profiles 
D(q). For statistical reasons leading to the above mentioned anomalies, we have resolved the D{q) profile only by 
three g-regions distinguishing between the three states: 1) mainly full-helix, 2) one and two-turn states, and 3) mainly 
unfolded states. We assume that these states are separated by the q- values, where the average helicity is 2/3 and 
1/3 in Fig. 2 (the salt-specific values of g2/3 and gi/3 being summarized in Tab. I). Although the resolution of D(q) 
is small we emphasize that it can almost quantitatively reproduce the folding and unfolding times plotted in Figs. 5 
and 6. Mean first passage times for folding and unfolding predicted by the FP equation involving salt-specific free 
energies and diffusivities are shown in the SI. 

Fig. 7 shows diffusivity profiles obtained by fitting to the average round-trip time defined in Eq. (0 (solid lines) 
and the Z)(g)-estimates resulting from fits to round-trip times of specific targets g t (symbols with error bars): first 
of all, we observe that they are not flat, a feature discussed previously for this peptide in the salt-free case£ The 
inhomogeneities reflect variations of the multi-dimensional configurational mobility of the peptide projected onto the 
one-dimensional RC q. Note also the fact that both analysis methods yield estimates, which coincide within error 
bars (the only exception being KC1 for q > gi/ 3 ), by this clearly validating our approach. After including the salt 
a few significant changes to D(q) are visible within the large error bars: firstly, in the solution with sodium salts 
there is a moderate decrease of the effective diffusivity by 30-60% in the large g-region, q > qi/3, where the peptide is 
mostly unfolded. In contrast, with KC1 the diffusivity seems to increase by 30-60% in the unfolded regions. Secondly, 
a drastic drop in diffusivity is observed for all salts in the central region g 2 /3 < g < gi/3, where the peptide features 
one or two a-helical turns; the decrease is about one order of magnitude for the sodium salts NaCl and Nal. Finally, 
we find a smaller decrease of the diffusivity again in the completely folded states q < g2/3, where the diffusion drops 
by about 30-60% for the KC1 and NaCl salts. 

Thus, while no clear-cut trends in the change of the diffusivity profile with salt can be recognized, clearly the 
diffusivity mainly decreases with salt, along with a significant drop in the partially folded states in the intermediate 
g-region. Importantly, the changes of D(q) with salt are obviously not just a rescaling of D(q) of the salt-free system 
as if the action would just stem from a nonspecific change due to a different bulk viscosity. In particular, the very 
few cases for which a viscosity argument applies could be for NaCl and KC1 in the low g-region, q < q^/3, or for NaCl 
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and Nal in the high g-region, q > qi/s' here the diffusivity is reduced by a factor of 2, roughly the same as for the 
viscosity increase (see Methods). Since there is no general trend, however, we must conclude that the change of the 
diffusivity profile originates from a combination of bulk viscosity effects and specific cosolute binding to the peptide. 
We believe that viscosity effects may be more important in the large q-region than in the low g-region: while large 
scale coil rearrangements in the solvent are frequent in the former, changes along the RC are mainly governed by 
internal mechanisms in the latter case; here only minor configurational rearrangements like the expulsion of one water 
molecule or one ion, or the forming and breaking of internal hydrogen-bonds take place. From this perspective it is 
interesting to see the diffusivity increasing in the high q-region for KC1, while from viscosity arguments only it should 
decrease by a factor of 2. 

However, given the diffusivity profiles in Fig. 7, the large increase of folding times [Figs 5 (a) and 6 (a)] for the 
sodium salts must thus be attributed not only to the changes in the free energy landscape but more dominantly, to 
the strongly reduced diffusivity in the intermediate q -region, q 2 /3 < q < Qi/3- The faster unfolding and the weaker 
dependence of unfolding vs. folding times on salt type (cf. Figs. 5 (b) and 6 (b)) seems to arise from a cancellation 
effect, where the small mobility in the intermediate q-regions is counterbalanced by the low unfolding barriers in F(q) 
(a comparison of the FP description for average folding and unfolding times with the raw data shown in Fig. 5 is 
given in the SI.) 

What are the molecular reasons for the major changes in the diffusivity profiles in the electrolytes solutions? Based 
on our structural and ion binding analysis in the previous section it is now easy to argue that the huge drop in the 
effective diffusivity in the sodium salts is generated by the long-lived configurations similar to those shown in Fig. 4. 
The long-lived character of these conformations is clearly observed in trajectory analysis and also manifested in the 
long binding times of cations on a nanosecond time scale shown in Fig. 3. While the form of the free energy landscape 
does not care about the life time of these states (just what fraction of time they are sampled), the long life times are 
clearly reflected in parts of the diffusivity profiles. Due to the stronger binding of Na + vs. K + to the peptide oxygen 
atoms the effect is much smaller in the KC1 solution than for the sodium salts. 



IV. SUMMARY AND CONCLUDING REMARKS 

In summary, we have investigated the specific effects of salt at molar concentration on the a-helical folding kinetics 
of a short, alanine-based and salt-bridge forming peptide by means of molecular simulations and a diffusional analysis. 
Mean folding times have been found to considerably depend on salt type with folding times varying over one order of 
magnitude. The molecular reason is the previously observed stronger binding affinity of Na + vs. K + ions to anionic 
peptide groups thereby transiently cross-linking multiple groups in the peptide. These binding processes increase 
the internal friction and induce a new, slow time scale. Within an analysis in terms of an effective diffusivity in a 
one-dimensional free energy landscape, these new time scales are expressed by a strong and salt-specific variation of 
the local diffusivity. A recent simulation study of a fully charged polyglutamic acid chain in salt solution showed that 
segmental relaxation kinetics were significantly slowed down due to the same molecular mechanisms^ 

The picture emerges that adsorption of ions not only alter the equilibrium but also kinetic properties of protein 
folding by direct binding mechanisms. Whether a general relation between preferential adsorption^, and changes in 
kinetics can be drawn may be an interesting notion for further research. Given the current insights it seems likely, 
however, that the change in kinetics not only depends on the amount of adsorbed ions but on the individual ion-peptide 
interactions. 

As we have demonstrated, molecular simulations can provide valuable information to understand the complex 
mechanisms in solvent-protein interactions and thereby protein stability and folding. The molecular mechanism 
found may be of general importance to understand cosolute effects on protein folding kinetics and shed more light onto 
experimentally observed cation-specific slowing down of (bio)poly-electrolyte kinetics ; 18 ' 21 in particular, for halophilic 
proteinS) 19 i 20 similar mechanisms may be at work in polymer melts^ More experimental studies are highly desirable; 
in particular, the novel long-lived loop-forming configurations in the denatured/unfolded states, in which sodium 
or similarly strong binders are bound and immobilized by the peptide backbone, may be experimentally accessible 
by nuclear magnetic relaxation dispersion methods (NMRDp or time-resolved FRET measurement^ probing salt- 
specific peptide relaxation and kinetics. 

Furthermore, the action of complex denaturants such as guanidinium and urea deserve further attention, and 
systematic studies on specific salt effects should follow. The guanidinium cation, for instance, has been shown to 
decrease friction in neutral (GlySer) n peptidesi 11 ' 12 We expect also a strong influence of other specifically binding 
cations on anionic peptides, such as lithium, or polyvalent cations, such as Mg +2 or Ca 2+ . Large effects may also be 
anticipated by exchanging the anion that has been found to considerably alter the unfolding kinetics of a halophilic 
protein i 19 ' 20 
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Fig. 1: Time series data for reaction coordinate q from the MD simulations in explicit salt water with varying salt 
type. Data in black shows the full 20 ps resolution, while data in red is smoothed over time windows of 2 ns. (figure 
not enclosed) 



Fig. 2: (a) Free energies F(q) for the EK peptide in different salt solutions; the profiles are shifted vertically 
for better comparison, (b) Average a-helicity of the peptide resolved by q. The snapshots illustrate the backbone 
structure of partially folded/unfolded states corresponding to values of the helicity indicated by arrows. Simulation 
snapshots are visualized using VMEfi. 



Fig: 3: Radial distribution function of cations around peptide oxygen atoms from the side chain (sc) carboxylates 
or backbone (bb) carbonyls. Inset: residence time distribution for the cations in the first solvation shell of the peptide 
oxygen atoms. The distribution for K + can be fitted by a single exponential with a time constant of r = 50 ps, while 
the one for Na + obeys a stretched exponential with r = 50 ps and stretching exponent j3 — 0.55. 



Fig. 4: Simulation snapshots of backbone configurations, in which a single Na + -ion (blue sphere) is trapped by 
multiple oxygen binding sites (red). These snapshots are dynamically selected such that the shown configurations 
existed for longer than 2 ns. In (a) and (b) no side chains are shown to better illustrate the binding to the backbone. 
In (c) and (d) all side chains are shown. Here, also glutamic acid side chains are involved in binding the cation, (figure 
not enclosed) 



Fig. 5: (a) Mean first passage times T{ for folding from state q — > q\. (b) Mean first passage times r u f for unfolding 
from state q — > q^. 



Fig. 6: (a) Mean first passage times T{ for folding to state q — > q%. (b) Mean first passage times r u f for unfolding to 
state g — >■ <72- 



Fig. 7: (a) Free energy landscapes (same as in Fig. 2 (a)), (b) Diffusivity profiles D(q) for all investigated systems 
resolved by three regions q < Q2/3> 92/3 < Q < Qi/Si an d q >i/3, which correspond to the fully helical, two-turn, 
and one-turn as well as unfolded states, respectively. Symbols including error bars result from fits to the ensemble 
of round-trip times TRT(q,qt) for different targets while fitting results to the average round-trip time curve trt 
(Eq. [7]) are displayed as solid lines. 



Fig. 8: Table of contents (TOC) figure, (figure not enclosed) 
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Figure 5 
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Figure 6 




